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BODIES  OF  REVOLUTION  HAVING  MINIMUM  DRAG 

AT  HIGH  SUPERSONIC  AIRSPEEDS1 

By  A.  J.  Eggers,  Jr.,  Meyer  M.  Resnikoff, 
and  David  H.  Dennis 


SUMMARY 


Approximate  shapes  of  nonlifting  bodies  having  minimum  pressure  fore¬ 
drag  at  high  supersonic  airspeeds  are  calculated.  With  the  aid  of  Newtons 
law  of  resistance,  the  investigation  is  carried  out  for  various  combina¬ 
tions  of  the  conditions  of  given  body  length,  base  diameter,  surface  area, 
and  volume.  In  general,  it  is  found  that  when  body  length  is  fixed,  the 
body  has  a  blunt  nose;  whereas,  when  the  length  is  not  fixed,  the  body 
has  a  sharp  nose.  The  additional  effect  of  curvature  of  the  flow  over 
the  surface  is  investigated  to  determine  its  influence  on  the  shapes  for 
minimum  drag.  The  effect  is  to  increase  the  bluntness  of  the  shapes  in 
the  region  of  the  nose  and  the  curvature  in  the  region  downstream  of  the 
nose.  These  shape  modifications  have,  according  to  calculation,  only  a 
slight  tendency  to  reduce  drag. 

Several  bodies  of  revolution  of  fineness  ratios  3  and  5,  including 
the  calculated  shapes  of  minimum  drag  for  given  length  and  base  diameter 
and  for  given  base  diameter  and  surface  area,  were  tested  at  Mach  numbers 
from  2.73  to  6.28.  A  comparison  of  theoretical  and  experimental  foredrag 
coefficients  indicates  that  the  calculated  minimum-drag  bodies  are  reason¬ 
able  approximations  to  the  correct  shapes.  It  is  verified,  for  example, 
that  tne  body  for  a  given  length  and  base  diameter  has  as  much  as  20  per¬ 
cent  less  foredrag  than  a  cone  of  the  same  fineness  ratio. 


INTRODUCTION 


The  shapes  of  nonlifting  bodies  of  revolution  having  minimum  pressure 
drag  at  supersonic  speeds  have  been  the  subject  of  numerous  theoretical 
investigations.  Karman  (reference  l)  determined  the  shape  of  such  a  body 
(neglecting  base  drag)  with  given  length  and  base  diameter.  Somewhat 

^-Supersedes  NACA  RM  A5IK27,  "Bodies  of  Revolution  for  Minimum  Drag  at 
High  Supersonic  Airspeeds,"  by  A.  J.  Eggers,  Jr.,  David  H.  Dennis, 
and  Meyer  M.  Resnikoff,  1952,  and  NACA  RM  A52D24,  "Supplementary  Note 
on  Modified- Impact-Theory  Calculations  for  Bodies  of  Revolution  Having 
Minimum  Drag  at  Hypersrnic  Speeds,"  by  Meyer  M.  Resnikoff,  1952. 
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later  Haack  (reference  2),  Ferrari  (reference  3),  Lighthill  (reference  4), 
and  Sears  (reference  5)  calculated  body  shapes  having  minimum  pressure 
drag  for  various  other  given  conditions  using  methods  similar  to  those 
first  employed  by  Karman.  In  all  these  investigations  the  assumption  of 
small  perturbation,  potential  flow  was  made.  It  is  to  be  expected,  there¬ 
fore,  that  the  shapes  obtained  by  these  investigators  are  representative 
of  minimum-drag  body  shapes  of  practical  fineness  ratios  at  low  supersonic 
Mach  numbers. 

Perhaps  the  first  calculation  of  the  shape  of  a  body  having  minimum 
drag  was  made  by  Newton  (reference  6)  using  a  method  analogous  to  the 
present  day  calculus  of  variations.  Newton  was  concerned  with  determining 
the  body  of  given  length  and  base  diameter  having  minimum  resistance  when 
moving  at  sufficiently  high  speeds  to  insure  that  the  inertia  forces  are 
large  compared  to  the  elastic  forces  in  the  immersing  fluid.  Thus,  as 
shown  by  Sanger  (reference  7)  and  Epstein  (reference  8),  the  law  of  resist¬ 
ance  adopted  by  Newton  approximates  that  (neglecting  viscous  forces)  for 
hypersonic  air  flows.  According  to  this  law,  the  local  resisting  pressure 
is  proportional  to  the  square  of  the  free- stream  velocity  component  normal 
to  the  body  surface.  Legerdre  (see,  e.g.,  reference  9)  further  investi¬ 
gated  Newton1  s  problem  and  concluded  that  if  no  restrictions  were  imposed 
on  the  variation  of  slope  along  the  surface,  a  body  having  a  meridian  curve 
composed  of  jagged  lines  (sharp  edges  forward)  could  be  constructed  which, 
with  this  law  of  resistance,  would  have  less  drag  than  Newton* s  body.  It 
may  easily  be  deduced,  however,  that  Newton’s  law  of  resistance  would  not 
be  satisfied  on  the  surface  of  Legendre’s  body  since  gas  would  be  trapped 
in  a  number  of  regions  along  the  jagged  contour.  It  may  be  shown  in  fact 
that  when  this  law  of  resistance  is  satisfied  at  the  surface  -  in  which 
case  the  surface  angles  must  lie  between  0  and  jt/2  radians  -  then 
Newton’s  body  may  be  considered  the  minimum  pressure  drag  body  for  the 
given  conditions . 

It  has  been  undertaken  in  the  present  report,  using  Newton’s  law  of 
resistance  and  the  calculus  of  variations,  to  determine  body  shapes  having 
minimum  pressure  drag  (neglecting  base  drag)  at  high  supersonic  speeds  for 
various  combinations  of  the  conditions  of  given  length,  base  diameter, 
surface  area,  and  volume.  The  effect  of  curvature  of  the  flow  over  the 
surface  is  also  investigated  to  determine  its  influence  on  the  shapes  for 
minimum  drag. 

Several  bodies  of  revolution,  including  two  of  the  bodies  determined 
from  this  analysis,  were  tested  at  Mach  numbers  from  2.73  t°  6.28  in  the 
Ames  10-  by  14-inch  supersonic  wind  tunnel .  Foredrag  data  at  zero  lift 
obtained  from  these  tests  are  compared  with  the  analytic  predictions  to 
assess  the  accuracy  of  the  theoretical  considerations. 
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local  cross-sectional  area  of  body 
local  speed  of  sound 
drag  coefficient 


(30 

(^) 


pressure  coefficient 

constant  of  integration 

pressure  foredrag 

maximum  body  diameter 

integrand  function 

drag  parameter  (  — P  -  ) 

\2tfqc  / 

hypersonic  similarity  parameter 


AO 


body  length 
Mach  number 


(?) 


distance  measured  normal  to  surface  of  body 

exponent  in  equation  defining  shapes  of  experimental  test  bodies 


static  pressure 
dynamic  pressure 

radius  of  curvature  of  streamline  in  plane  containing  axis  of 
symmetry  (i.e.,  meridian  plane)  of  body 

body  surface  area 

resultant  Telocity 

body  volume 

coordinates  of  point  on  meridian  curve  of  body  (origin  of  coordinate 
system  coincides  vith  nose  of  body,  and  X  axis  coincides  with 
axis  of  symmetry) 


4 


NA.CA.  TN  3666 


7  ratio  of  specific  heat  at  constant  pressure  to  specific  heat  at 
constant  Yolume 

6  angle  (in  meridian  plane)  between  free- stream  direction  and 
tangent  to  body  surface 

X  Lagrange  multiplier 

p  density 

Subscripts 

o  free-stream  conditions 

1  values  at  nose  point  of  meridian  curve 

2  values  at  base  point  of  minimizing  curve 

+  right-hand  limiting  value  of  quantity  at  corner  on  minimizing  curve 
—  left-hand  limiting  value  of  quantity  at  comer  on  minimizing  curve 
B  values  along  meridian  curve 
c  cone  values 


THEORY 


The  investigation  undertaken  here  is  concerned  with  the  shapes  of 
nonlifting  bodies  of  revolution  having  mini  mum  pressure  foredrag  at 
high  supersonic  airspeeds.  Difficulties  inherent  in  the  calculation  of 
these  shapes  make  it  desirable  to  simplify  the  drag  equation  insofar  as 
is  practicable,  consistent  with  retaining  the  salient  features  of  the 
dependence  of  drag  an  body  shape  and  free-stream  conditions.  Likewise, 
in  view  of  the  several  conditions  to  be  treated  (viz.,  given  length, 
base  diameter,  volume,  and  surface  area),  it  is  convenient  to  set  up  a 
procedure  of  analysis  to  fit  the  general  problem  at  hand.  These  funda¬ 
mental  considerations  will  be  discussed  prior  to  the  determination  of 
specific  minimum-drag  shapes. 
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Fundamental  Considerations 


Simplified  drag  theory.-  As  pointed  out  in  the  introduction, 
Newton’s  lav  of  resistance  applies  approximately  to  bodies  traveling 
at  high  supersonic  airspeeds.  This  observation  has  basis  in  the  fact 
that  at  such  speeds  the  inertia  forces  predominate  over  the  elastic 
forces  in  the  disturbed  air.  Thus,  oblique  shock  flows  approach  the 
corpuscular- type  flows  treated  by  Newton  as  the  Mach  number  of  the  free 
stream  becomes  large  compared  to  1.  If  it  is  further  assumed  that 
7  of  the  disturbed  fluid  approaches  l,2  the  shock-wave  angle  approaches 
the  flow-deflection  angle  (see  sketch).  In  this  case  the  pressure 


U 


Q. 


coefficient  at  a  point  Just  downstream  of  the  wave  is  given  by  the 
simple  expression  (reference  8) 


Cp  =  2  sin2  5  (l) 

This  equation  is  recognized,  of  course,  as  being  (aside  from  the  con¬ 
stant  multiplier)  a  mathematical  statement  of  Newton’s  law  of  resist¬ 
ance  for  corpuscular  or  impact— type  flow. 

When  the  curvature  of  the  body,  and  hence  of  the  disturbed  flow, 
is  small  in  the  stream  direction,  equation  (l)  should  also  predict  the 
pressure  coefficients  at  the  surface  of  a  body  since,  in  this  case,  the 
centrifugal  forces  in  the  thin  layer  of  air  (sometimes  referred  to  as 
the  hypersonic  boundary  layer)  between  the  shock  and  the  surface  should 
not  appreciably  alter  the  impact  pressures.  When  the  curvature  of  the 
body  is  large  in  the  stream  direction,  centrifugal  forces  in  the  fluid 
between  the  shock  «md  the  surface  may  alter  the  pressures  at  the  surface 

2This  assumption  was  perhaps  first  suggested  by  Epstein  (reference  8). 

It  has  the  advantage  of  simplifying  the  impact  force  analysis  without 
significantly  influencing  its  accuracy.  In  other  words,  these  forces 
in  flows  about  axially  symmetric  bodies  are  relatively  insensitive  to 
changes  in  7  from  1.4  to  1.  Interestingly  enough,  present  indica¬ 
tions  are  that  disturbed  air  flows  in  flight  at  extremely  high  Mach 
numbers  will  be  characterized  by  values  of  7  below  1.4  (see  ref¬ 
erences  10  and  11). 
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from  those  Just  downstream  of  the  shock,  Busemann  (reference  12) 
investigated  this  problem  and  found  that  the  pressure  coefficient  at 
a  point  on  the  surface  of  a  body  curved  in  the  stream  direction  is 
given  by  the  relation 

Cp  =  2^sin2  5  +  sin  5  ~  J'  cos  5  dA  ^  (2) 

'  o 

in  the  limit  as  M— >  »and  7  — >  1. 

In  order  to  assess  the  accuracy  with  which  the  preceding  equations 
may  be  expected  to  provide  the  pressure  distributions,  and  thus  pres¬ 
sure  drags,  on  bodies  operating  at  high  supersonic  airspeeds,  the  pre¬ 
dictions  of  these  equations  are  compared  in  figure  1  with  those  of  the 
method  of  characteristics  (obtained  from  reference  13)  for  an  ogive 
operating  at  a  value  of  the  hypersonic  similarity  parameter  K  (ratio 
of  free-6tream  Mach  number  to  slenderness  ratio)  equal  to  2,  corre¬ 
sponding  to  a  free— stream  Mach  number  of  6.  It  is  evident  that  the 
theory  of  Busemann  (equation  (2))  yields  far  too  low  pressures  down¬ 
stream  of  the  nose,  while  the  simple  impact  theory  (equation  (l))  is  in 
reasonably  good  over-all  agreement  with  the  method  of  characteristics.3 
The  relatively  poor  predictions  of  the  Busemann  theory  are  associated 
with  the  fact  that  it  strongly  overestimates  centrifugal— force  effects 
at  free-6tream  Mach  numbers  which  are  large  compared  to  1,  but  for  which 
7  of  the  air  flow  downstream  of  the  bow  shock  is  closer  to  1.4  than  1 
(i.e.,  at  flow  conditions  of  principal  interest  in  this  paper).  This 
nuitter  will  be  discussed  in  greater  detail  later  in  the  paper.  Agree¬ 
ment  comparable  to  that  Just  discussed  is  obtained  with  the  other 
results  resented  in  reference  13  for  K  =  2.  For  lower  values  of  K 
the  agreement  of  the  impact  theory  with  the  method  of  characteristics  is 
somewhat  poorer,  as  would  be  expected;  however,  it  does  not  become 
unacceptably  poor  except  for  values  of  K  below  1  (e.g. ,  the  pressure 
coefficients  differ  by  from  0  to  35  percent  for  a  K  of  l/2).  It  is 
therefore  concluded  that  for  values  of  K  greater  than  1,  equation  (l) 
may  be  used  with  acceptable  accuracy  for  the  purposes  of  this  paper  to 
predict  the  pressure  distributions  and  thus  pressure  drags  on  bodies.4 

3The  characteristics  solutions  of  reference  13  were  carried  out  for 
7  »  1.4.  Free  flight  at  the  larger  values  of  K  considered  in  this 
report  would  produce  values  of  7  downstream  of  the  bow  shock  slightly 
less  than  1.4;  however,  the  small  decreases  in  7  would  not  signifi¬ 
cantly  alter  the  results  presented. 

'‘The  approximate  theory  of  reference  14  might  also  be  employed  in  an 
analysis  of  this  type  —  this  theory  is,  in  fact,  more  accurate  than 
those  discussed,  so  long  as  the  flow  is  everywhere  supersonic.  How¬ 
ever,  the  theories  considered  here  have  the  advantage,  as  will  be 
apparent  later,  of  predicting  approximately  correct  values  of  surface 
pressures  for  arbitrarily  large  surface  slopes. 
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For  this  reason,  and  "because  of  its  simplicity,  it  is  employed  through¬ 
out  the  subsequent  analysis. 

If  the  manner  in  which  the  pressure  coefficient  varies  over  the 
surface  is  known,  it  is  a  simple  matter,  of  course,  to  evaluate  the 
pressure  drag  of  a  bu  .j„  Neglecting  the  base-drag  contribution,  we 
have  then 


(L.  q  rtd2  n 

D  *  -2 — 2 - *  2*  qQ  /  Cpyy*dx 

where  y*  denotes  the  derivative  dy/dx.  This  equation  may  be 
expressed  in  a  form  more  convenient  for  use  here 


(3) 


w 


If  Gp  in  this  expression  is  replaced  by  its  value  given  in  equation  (l) 

r»2 


(noting  that  sin2  5  « 


_  y1 


1+y* 

1^ 


-g),  there  i6  then  obtained  the  relation 

J. 


2yy* 


dx 


1+y* 


(5) 


It  remains  now  to  consider  the  procedure  for  employing  this  expression 
in  combination  with  the  methods  of  the  variational  calculus  in  order  to 
determine  the  desired  minimum-drag  body  shapes. 


Procedure  for  Calculating  Minimum-Drag  Bodies 


The  calculation  of  minimum-drag  body  shapes  of  interest  here  is 
equivalent  to  determining  the  form  of  the  function  y  *  y(x)  which 
minimizes  the  integral  defined  in  equation  (5)  for  the  various  given 
conditions.  In  considering  the  procedure  for  carrying  out  this  calcu¬ 
lation,  however,  it  is  convenient,  for  reasons  that  will  be  apparent 
later,  to  write  equation  (5)  in  a  form  which  effectively  yields  the 
total  drag  as  the  sum  of  the  drag  on  any  finite  region  of  infinite  slope 
at  the  nose  plus  the  drag  on  the  surface  downstream  of  the  nose.  Thus 
we  have 


Id  =  7i2  + 


guilds 

1+y*2 


(6) 
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where  the  variable  limit  x2  is  introduced  to  permit  variations  in 
body  length.  The  conditions  of  given  volume  or  given  surface  area  are 
fixed  hy  the  auxiliary  requirements  that,  respectively, 

x2 

X  r  y2  dx  *  const.  (7) 

* 

or  (neglecting  base  area) 

yi2  p 

~  * - +  /  y */l  +  y,s  dx  =  const.  (8) 

2iC  2  J 
o 

When  the  length  and  base  diameter  are  given,  the  problem  is  simply 
to  minimize  the  function  ID  given  by  equation  (6).  However,  according 
to  the  isoperimetric  rule  of  the  calculus  of  variations  (see,  e.g. , 
reference  15),  the  problem  of  minimizing  the  function  subject  to 

the  auxiliary  condition  given  by  equation  (7)  or  (8),  is  equivalent  to 
minimizing  the  new  function  JD,  where 

JD  =  ^  +  X  7  (9> 

or 

jd  -  xd  +  x  h  (10) 

depending  on  whether  the  volume  or  surface  area  is  given.  The  parameter 
X  is  a  constant,  sometimes  called  the  Lagrange  multiplier. 

With  the  aid  of  equations  (6)  through  (10),  the  integrand  functions 
to  be  minimized  can  be  immediately  written.  These  functions  are  as 
follows : 


case  a,  given  length  and  base  diameter 

_  2yy* 


i+y 


T5" 


case  b,  given  volume  and  length  or  base  diameter 

f  =  Sail  +  Xy2 
1+y'2 


(11) 


(12) 
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case  c,  given  surface  area  and  length  or  base  diameter 

f  *  ^L'--  +  Xy^/l  +  y'2  (13) 

1+y* 

How  any  function  y  *  y(x)  which  minimizes  equation  (6),  (9),  or 
(10)  must,  irrespective  of  the  given  conditions,  satisfy  the  Euler 
equation  (for  zero  first  variation  of  ID  or  Jjj  with  small  changes  in 
the  function  y(x)) 

5  V  -  fy  "  0 

where  f  '  and  f„  denote  the  partial  derivatives  and  respec— 

*  *  a  7  dy 

tively.  Since  the  integrand  functions  given  above  are  free  of  the 
independent  variable,  the  first  integral  of  the  Euler  equation  for 
these  functions  follows  immediately,  namely, 

y1  fy  1  —  f  *  const.  (15) 


Substituting,  successively,  equations  (ll),  (12),  and  (13)  into  this 
equation  there  are  then  obtained  the  expressions 


and 


%y' 


(l+y,2)‘ 


const. 


4yy» 


(i+y,2)‘ 


—  Xy2  =  const. 


4y* 


L(i+y'2)2  */i+y,5j 


const  • 


(16) 

(17) 

(18) 


for  cases  a,  b,  and  c,  respectively.  Solutions  to  these  differential 
equations  satisfying  the  terminal  conditions  on  the  bodies  are  minimiz¬ 
ing  curves  for  the  given  conditions. 

When  the  end  points  of  a  minimizing  curve  are  not  fixed,  other 
terminal  conditions  must  be  imposed  on  the  function  y  *  y(x) .  Thus, 
to  determine  the  ordinate  at  the  nose,  it  is  required  that  (see  ref¬ 
erence  16) 
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f  -  JL  y2 

dy  ^ 


for  cases  a  and  b,  while 


7=71 


gy  y’2g=# 
(i+y'  ) 


=  0 


(19) 


y=yi 


L  V  ~  x  f  > 


7=*7i 


=  7l  [2  y'4+3y-'4  +  »^£L-(2+x)  =0 

L  (l+y*2)  Vl+7*  Jv=!V’ 


^=7i 


(20) 


for  case  c.  Similarly,  when  the  length  is  not  given  it  is  necessary 
that 


y'fy.  -  f 


=  0 


(21; 


x=x2 


and  when  the  "base  diameter  is  not  given  it  is  required  that 


V 


=  0 


|X=X2 


(22) 


In  addition  to  the  above  described  conditions,  two  checks  must  be 
made  to  determine  completely  the  shape  of  a  minimizing  curve.  The  first 
of  these  checks  entails  ascertaining  whether  there  are  any  corners 
(between  the  end  points)  on  the  curve.  This  is  accomplished  by  deter¬ 
mining  whether  the  function  y  =  y(x)  can  satisfy  the  requirement  that 
(see  reference  15) 


(23) 


at  a  point  of  discontinuity  in  y'.  If  this  equation  is  not  satisfied, 
no  comers  exist.  The  second  check  requires  that  the  Legendre  condition 
(for  a  positive  second  variation), 

*rr  >  0  (2*0 

be  satisfied  everywhere  on  the  curve.  With  the  aid  of  these  checks,  the 
minimizing  curves  far  various  combinations  of  the  conditions  of  given 
length,  base  diameter,  volume,  and  surface  area  can  be  uniquely  defined. 
The  calculation  of  these  curves  for  several  such  combinations  is  now 
undertaken. 
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Calculation  of  Minimum-Drag  Bodies 


Given  length  and  base  diameter.-  Equations  (l6)  and  (19)  give  the 
first  integral  to  Euler’s  equation  and  the  terminal  condition  at  the 
nose,  respectively,  for  these  given  conditions.  It  is  evident  upon 
examining  these  equations  that  the  minimizing  curve  cannot,  in  general, 
pass  through  both  the  points  (0,0)  and  (x2,y2),  hut  must,  in  fact,  have 
its  forward  termination  point  at  (o,yi)  withy!*  *  1.  With  this  infor¬ 
mation,  the  minimizing  curve  can  he  represented  in  parametric  form, 
namely. 


_  yi  (1  +  y,2)‘ 


-2—  +  -  I  +  In  y' 


It  is  easily  shown  with  the  solution  to  the  Euler  equation  and  equar- 
tion  (23)  that  there  are  no  corners  on  the  minimizing  curve;5  thus  the 
variation  of  y  with  x  is  readily  determined  with  the  relations  of 
equation  (25)  for  a  given  l  and  d  (corresponding  to  a  given  x2  and  y2) 
of  a  hody.  These  relations  for  a  body  of  given  fineness  ratio  can  he 
shown  to  he  equivalent  to  those  originally  developed  hy  Newton  (see  ref¬ 
erence  6). 

Given  length  and  volume.—  For  these  given  conditions,  the  terminal 
conditions  (equations  (19)  and  (22))  require  the  slopes  at  the  nose  and 
at  the  base  to  he,  respectively,  yi*  *  1  and  y2*  =  0.  The  first  inte¬ 
gral  to  the  Euler  expression  (equation  (17))  then  leads  to  the  following 
parametric  representation  of  the  minimizing  curve: 


X(l+y'2)” 


r7  dy 

I=J  r 

7i  J 


gy* 

-x(l+y>2) 


yi  -  Xyi 


Similarly,  it  can  he  shown  that  there  are  no  corners  between  (o,yi] 
and  (x2,y2)  on  any  of  the  minimizing  curves  to  he  treated  here. 
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From  the  relations  of  equation  (26)  it  is  clear,  again,  that  the  mini¬ 
mizing  curve  cannot  pass  through  (0,0),  the  condition  yi1  =  1  deter¬ 
mining  a  value  yi  >0.  These  relations,  together  with  the  volume 
condition  (equation  (7))  and  the  given  length  condition,  serve  to 
determine  yi  and  X  and  thus,  of  course,  the  shape  of  the  entire  body. 
As  the  length  approaches  0,  X  becomes  infinitely  negative;  vhile,  as 
the  length  becomes  infinitely  large,  X  approaches  0.  (in  the  latter 
case  the  body  shape  approaches  the  minimum-drag  shape  for  the  given 
length  and  diameter  condition,  l/d  -$>*.)  Intermediate  negative  values 
of  X  correspond  to  intermediate  values  of  length  for  a  given  volume. 


Given  length  and  surface  area.—  In  this  case  a  first  integral  to 
the  Euler  equation  is  given  by  equation  (l8),  and  the  parametric  repre¬ 
sentation  of  the  minimizing  curve  may  be  -written  immediately  in  the 
form 


_  const,  (l+y*2) 

7  3  2  3/2 

4y»  —  X(l+y*  ) 

Z-  f  % 

J7i  y' 


(27) 


Upon  examination  of  this  equation  and  equations  (20)  and  (22),  it 
becomes  apparent  that,  again,  the  minimizing  curve  cannot  go  through 
the  point  (0,0).  The  latter  equations  determine  uniquely,  however,  the 
values  of  yi*  (yi*<  1),  and  y2»  (o<y2,<yi,)>  in  terms  of  the  param¬ 
eter  X.  Similarly,  the  length  and  surface-area  condition  in  combinar- 
tion  with  the  above  equations  determines  the  value  of  X.  Thus  it  is 
easily  shown  that  the  practical  range  of  X  is  from  -2  to  0  (corre¬ 
sponding  to  body  lengths  of  from  zero  to  infinity  for  a  given  surface 
area  —  in  the  latter  case  the  Newton  body  is  again  obtained). 


Given  base  diameter  and  volume.—  With  these  given  conditions,  the 
first  integral  to  the  Euler  relation  is  given  by  equation  (17),  while 
the  terminal  conditions  at  the  fore-and-aft  ends  of  the  body  are  fixed 
by  equations  (19)  and  (21),  respectively.  It  is  evident  that  the  mini¬ 
mizing  curve  must,  in  general,  pass  through  the  origin  in  order  to 
satisfy  all  these  equations  in  addition  to  the  Legendre  condition 
(equation  (24)).  The  shape  of  the  minimizing  curve  may  thus  be  defined 
parametrically  as  follows: 
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7 


4  7'3 

x  (l+y*2)2 


x  g  y'4  +  3r,a 
^  (l+y'2)2 


> 


(28) 


where  y^ 1  =0.  Combining  this  expression  with  equation  (7)>  there  is 
then  obtained  for  the  volume  of  the  body 


V  =  JIzL-  (y2t4  +  6y2'2  +  45)  (29) 

120  y2» 

The  range  of  X  for  vhich  these  results  are  applicable6  is  from 

zero  to  3</3/4y2>  corresponding  to  a  volume  range  from  infinity  to 

it  y2SvT/5.  For  a  given  y2  and  a  given  V>*  y23  /J/5  (corresponding 

to  ^>*/J/2),  equation  (29)  has  tvo  solutions  in  y2*.  One  solution 

yields  values  of  y2*  greater  than  «/3,  which  result  violates  the 
Legendre  condition  (see  equation  (24)),  while  the  other  yields  permis¬ 
sible  values  less  than  </3.  When  y2  and  y2 ’  are  known,  X  may  then  be 
determined  from  the  first  relation  of  equation  (28),  namely, 


X 


—  ya'3 

ya  (l+y2'2)5 


(30) 


The  determination  of  y  and  x  follows  directly,  of  course,  from 
equation  (28). 

Given  base  diameter  and  surface  area.—  In  this  case  equations  (l8), 
(20),  and  (21)  determine  the  shape  of  the  minimizing  curve  as  being 
simply  a  straight  line 

y  =  x 


(4/4) 


2/3 


1  -  (4/4) 


2/3 


(31) 


where  the  parameter  X  is  given  by  the  equation 

6Ibe  solution  given  here  is  not  applicable  to  bodies  of  extremely  small 
fineness  ratios  (viz,,  ^<^)  as  can  be  easily  deduced  from 
equation  (28). 
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X  -  h  (ity22/s)3  (32) 

Thcs,  the  minimum- drag  body  for  given  "base  diameter  and  surface  area  is 
a  con?. 


Comparison  of  Minimum-Drag  Body  Shapes 


The  previous  calculation  of  minimum-drag  bodies  reveals  two 
general  characteristics  of  their  shapes;  namely,  when  the  length  is 
given  (fixed)  the  bodies  assume  blunt  noses,  whereas,  when  the  length 
is  not  given  (i.e.,  is  free),  the  bodies  assume  sharp  noses.  The 
former  characteristic  may  be  traced  to  the  fact  that  with  the  length 
restricted,  the  net  drag  is  reduced  by  accepting  higher  pressures  on  a 
relatively  small  area  of  large  slope  near  the  nose,  thus  achieving 
lower  pressures  on  a  relatively  large  area  of  small  slope  near  the  base. 
On  the  other  hand,  when  the  length  is  not  restricted  it  is  evident  that 
a  sharp  rather  than  a  blunt  nose  will  obtain  for  minimum  drag,  since 
the  drag  of  any  blunt-nosed  body  can  be  reduced  by  simply  relaxing  the 
requirement  on  length,  thereby  allowing  the  body  to  be  made  sharp  nosed 
and  generally  more  slender. 

In  order  to  permit  a  quantitative  comparison  of  the  shapes  of  the 
calculated  minimum- drag  bodies,  typical  meridian  curves  for  these  bodies 
are  shown  in  figure  2.  For  simplicity  the  bodies  are  compared  on  the 
basis  of  the  same  fineness  ratio  —  ordinates  have  been  plotted  to  an 
expanded  scale  to  better  indicate  the  relative  shapes.  The  maximum 
bluntness  is  evidently  obtained  when  the  drag  is  minimized  for  a  given 
length  and  surface  area,  while  the  maximum  sharpness  (a  cusp  nose)  is 
obtained  when  the  base  diameter  and  volume  are  given.  It  is  apparent 
from  figure  2  that  the  flat-nosed  portions  of  the  meridian  curves  for 
the  given  length  bodies  are  in  all  cases  very  small.  For  example, 
yi  equals  0.0050ye  for  the  body  of  given  length  and  volume.  On  the  basis 
of  several  calculations  it  is  indicated,  as  might  be  expected,  that  the 
degree  of  bluntness  will  increase  with  decreasing  fineness  ratio. 

It  is  also  of  interest  to  compare  minimum-drag  body  shapes  deter¬ 
mined  with  the  aid  of  the  linear  theory  (see,  e.g. ,  reference  2)  with 
those  found  using  the  impact  theory,  that  is,  bodies  especially  suited 
for  flight  at  low  and  high  supersonic  speeds,  respectively.  Such  a 
comparison  is  shown  in  figure  3  for  the  case  of  given  length  and  base 
diameter.  It  is  seen  that  qualitatively  the  shapes  are  similar  although 
the  minimum-drag  body  for  low  supersonic  speeds  is  generally  the  fatter 
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of  the  two.7  Comparisons  of  the  results  of  this  paper  with  those  of 
reference  2  for  other  given  conditions  also  indicate  qualitative  agree¬ 
ment  as  to  general  body  shapes  despite  the  marked  difference  in  the  laws 
governing  the  surface  pressures. 

All  the  preceding  analysis  has  been  predicated  on  the  assumption 
that  the  flow  of  air  at  high  supersonic  speeds  may,  insofar  as  pressure 
forces  are  concerned,  be  approximated  by  a  Newtonian- type  flow.  It 
remains  now  to  test  the  accuracy  of  this  assumption  and  other  aspects 
of  the  analysis  by  experiment. 


EXPERIMENT 


It  has  been  undertaken  to  obtain  a  partial  check  on  the  findings  of 
the  preceding  theoretical  analysis  by  determining  experimentally  the  fore¬ 
drags  on  a  family  of  bodies  of  given  fineness  ratios  at  Mach  numbers  from 
2.73  to  6.28.  The  analysis  may  be  expected  to  apply,  at  least  approxi¬ 
mately,  in  this  range  since  for  the  bodies  tested  the  corresponding  values 
of  the  hypersonic  similarity  parameter  K  were,  for  the  most  part,  greater 
than  1.  A  brief  description  of  these  tests  is  now  presented. 


Apparatus  and  Tests 


The  tests  were  conducted  in  the  Ames  10-  by  1 4-inch  supersonic  wind 
tunnel,  which  is  of  the  continuous -flow  nonreturn  type  and  operates  with  a 
nominal  supply  pressure  of  6  atmospheres .  The  Mach  number  in  the  test 
section  may  be  varied  from  approximately  2.7  to  6.3  by  changing  the  rela¬ 
tive  positions  of  the  symmetrical  top  and  bottom  walls  of  the  wind  tunnel. 
During  operation  at  the  higher  Mach  numbers,  the  supply  air  is  heated 
before  it  enters  the  wind  tunnel  to  prevent  condensation  of  the  air.  A 
detailed  description  of  the  wind  tunnel  and  its  associated  equipment  and 
of  the  characteristics  of  the  flow  in  the  test  section  may  be  found  in 
reference  17. 

Aerodynamic  drag  forces  were  measured  with  a  strain-gage  balance. 

Tare  forces  on  the  sting  supports  were  essentially  eliminated  by  shrouds 
that  extended  to  within  0.040  inch  of  the  model  base.  Axial  forces  on  the 
bases  of  the  models  were  determined  from  measured  base  pressures  and  from 

7Part  of  this  difference  in  shapes  stems  from  the  fact  that  the  body 

derived  using  linear  theory  was  required  to  have  zero  slope  at  the  base. 
Also,  as  will  be  shown  later,  the  true  minimum-drag  shape  at  high  super¬ 
sonic  airspeeds  may  be  somewhat  fatter  than  that  obtained  using  impact 
theory,  due  to  the  fact  that  centrifugal  forces  are  neglected  in  this 
theory. 
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free-stream  static  pressures  and  were  subtracted  from  measured  total  drag 
forces;  thus,  the  data  presented  do  not  include  the  forces  acting  on  the 
bases  of  the  test  bodies. 

Reynolds  numbers  based  on  the  maximum  diameter  of  the  test  bodies 

were: 


Mach  number 

Reynolds  number 
million 

2.75 

0.70 

3.50 

.95 

4.00 

.72 

5.05 

•  35 

6.28 

.15 

Reynolds  numbers  based  on  model  length  may  be  obtained  by  multiplying  the 
above  values  by  model  fineness  ratio. 


Models 


Five  models  of  fineness  ratio  3,  and  three  models  of  fineness  ratio  5 
were  tested.  With  the  exception  of  an  z/d  =  3  tangent  ogive  (this  shape 
was  included  as  being  typical  of  those  in  common  usage),  all  models  had 
meridian  section  shapes  given  by  the  equation 

*  ■  (f)"  <»> 

where  n  was  given  values  of  1,  3/4,  1/2,  and  l/4.  When  n  =  3/4,  the 
body  shapes  defined  by  the  above  expression  closely  approximate  the 
minimum-drag  shapes  for  given  length  and  base  diameter  (equation  (25) )8 
for  Z/d  =  3  and  5  (see  fig.  4). 

When  n  =  1,  the  cone  is,  of  course,  obtained  which  is  the  minimum- 
drag  body  for  a  given  base  diameter  and  surface  area.  Minimum-drag  shapes 
for  two  different  given  conditions  are  thus  included  among  the  bodies 
tested. 

Photographs  of  the  eight  models  tested  are  shown  in  figure  5.  The 
Z/d  =  3  bodies  (fig.  5(a))  are,  from  left  to  right  in  the  photograph, 
the  cone,  3/4-power  body,  l/2-power  (parabolic)  body,  l/4-power  body, 
and  the  tangent  ogive  \*hich  has  a  profile  section  radius  of  curvature 
of  9*25  body  diameters.  From  left  to  right  in  figure  5(b)  are  the 
Z/d  =  5  cone,  3/4-power  body,  and  l/2-power  body.  The  base  diameter  of 
all  models  was  1  inch. 

“The  accuracy  of  this  approximation  increases  with  increasing  values  of 
Z/d  as  can  easily  be  seei.  upon  examination  of  equation  (25) » _ 
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Accuracy  of  Test  Results 

The  accuracy  of  the  foredrag  coefficients  is  effected  by  uncertain¬ 
ties  in  the  measurements  of  the  following  quantities:  stagnation  pres¬ 
sures,  free-stream  static  pressures,  base  pressures,  and  the  forces  on 
the  models  as  measured  by  the  strain-gage  balance. 

Both  static  and  free-stream  dynamic  pressures  were  determined  from 
wind-tunnel  calibration  data  and  stagnation-pressure  readings.  The  latter 
measurements  were  accurate  to  within  ±l/2  percent,  thus  static  and  dynamic 
pressures  are  uncertain  by  this  amount  plus  possible  calibration  errors 
of  ±1  percent  over  the  Mach  number  range  of  the  tests .  The  uncertainty 
in  foredrags  due  to  inaccuracies  in  the  determination  of  base  pressures 
does  not  exceed  ±1  percent. 

Because  of  the  small  drag  forces  measured,  the  source  of  greatest 
error  was  the  strain- gage  balance  system.  The  uncertainty  in  drag  due 
to  zero  shifts,  thermal  effects,  and  friction  varied  from  approximately 
±2  percent  at  the  lower  Mach  numbers  to  ±6  percent  at  the  highest  Mach 
numbers. 

The  combined  effects  of  all  the  sources  of  error  result  in  probable 
uncertainties  in  measured  foredrag  coefficients  of  from  ±0.001  at  the  low 
Mach  numbers  to  ±0.005  at  a  Mach  number  of  6.28.  In  order  to  reduce  this 
error  in  the  data  presented  here,  particularly  at  the  higher  Mach  numbers, 
several  measurements  were  made  at  each  Mach  number  and  the  average  values 
of  foredrag  coefficients  were  employed. 

RESULTS  AND  DISCUSSION 


The  variations  with  Mach  number  of  the  measured  foredrag  coefficients 
are  shown  in  figure  6.  It  is  evident  that  the  3/4-power  bodies  do,  as 
predicted,  have  the  minimum  foredrags  of  all  the  test  bodies  with  the  same 
fineness  ratio,  the  drag  of  the  3/4-power  body  being  as  much  as  20  percent 
less  than  that  of  the  cone  of  the  same  fineness  ratio.  The  general 
increase  in  foredrag  at  Mach  numbers  in  the  neighborhood  of  5  and  greater 
can  be  traced  to  an  increase  in  friction  drag.  This  latter  increase  is, 
in  turn,  caused  by  the  relatively  large  decrease  in  Reynolds  number  with 
increasing  Mach  number  in  this  range  (see  section  on  Apparatus  and  Tests). 

A  check  on  the  over-all  accuracy  with  which  the  optimum  shapes  are 
predicted  by  the  analysis  is  obtained  by  comparing  theoretical  and  experi¬ 
mental  values  of  the  relative  foredrag  coefficients  of  the  test  bodies . 

Such  a  comparison  is  given  in  figure  7  where  the  ratios  of  the  foredrag 
coefficients  of  a  test  body  to  the  corresponding  coefficients  of  the  cone 
of  the  same  fineness  ratio  are  shown  as  a  function  of  the  exponent  n  in 
equation  (33)  which  defines  the  shapes  of  the  test  bodies.  The  theoretical 
predictions  of  the  impact  theory  appear  to  be  in  good  agreement  with  the 
experimental  results  at  the  higher  values  of  n  (approximately  n  >  0.6). 
Thus  it  is  suggested  that  the  3/4-power  body  is  a  reasonable  approximation 
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to  the  correct  minimum- foredrag  shape  of  given  fineness  ratio.  At  the 
lower  values  of  n,  however,  it  is  indicated  that  the  relative  drag  is 
significantly  overestimated  by  this  theory.  This  result  is  not  entirely 
surprising  since  the  theory  neglects  centrifugal-force  effects  in  the 
disturbed  flow,  and  these  effects  must  appreciably  alter  the  pressures 
over  the  highly  curved  noses  of  the  blunter  bodies. 

As  discussed  earlier,  the  Busemann  theory  for  infinitely  high  Mach 
numbers  overestimates  these  effects  at  the  Mach  numbers  of  interest  here. 

It  has  therefore  been  undertaken  in  Appendix  A  of  this  paper  to  obtain  a 
better  estimate  of  centrifugal  forces  by  accounting  approximately  for  the 
decrease  in  these  forces  (at  finite  but  high  Mach  numbers)  associated  with 
the  increase  in  the  lateral  extent  of  the  distrubed  flow  field  with 
increasing  distance  downstream  from  the  nose  of  the  body.  The  predictions 
of  the  modified  impact  theory  shown  in  figure  7  were  obtained  with  the 
aid  of  this  estimated  centrifugal-force  effect  (see  equation  (A9))  in 
combination  with  equations  (l)  and  (3).  It  is  indicated  that  this  theory 
is  markedly  superior  to  the  impact  theory  at  the  lower  values  of  r.;  cor¬ 
responding  to  the  blunter  bodies,  over  the  test  Mach  number  range.  The 
estimate  of  the  centrifugal  forces  would  thus  appear  to  be  in  fair  agree¬ 
ment  with  the  actual  magnitude  of  these  forces. 

It  is  also  indicated  in  figure  7  that  at  the  higher  test  Mach  numbers 
(hence  higher  values  of  K)  ,9  the  modified  theory  is  generally  somewhat 
superior  to  the  impact  theory.  This  result  suggests  that  improved  approxi¬ 
mations  to  the  correct  minimum- foredrag  shapes  for  values  of  K  appreci¬ 
ably  greater  than  1  may  be  obtained  by  using  this  theory  rather  than  the 
simple  impact  theory.  Accordingly,  calculations  of  minimum-drag  shapes 
have  been  made  using  the  modified  impact  theory  in  the  manner  discussed 
in  Appendix  B.  The  body  shapes  obtained  (see  Appendix  B)  are  for  the  same 
given  geometric  conditions  as  those  previously  determined  using  impact 
theory.  The  resulting  shape  for  given  length  and  diameter  is  shown  in 
figure  8.  Newton* s  body  of  the  same  fineness  ratio  is  also  shown  for 
comparison.  The  body  shape  determined  by  the  modified  theory  is  somewhat 
more  blunt  in  the  region  of  the  nose  and  has  more  curvature  in  the  region 
downstream  of  the  nose  than  Newton* s  body.  A  similar  comparison  is  shown 
in  figure  9  for  the  bodies  of  given  base  diameter  and  surface  area.  In 
this  case  both  bodies  have  pointed  noses  because  the  length  is  not  fixed, 
but,  in  the  same  manner  as  for  the  bodies  of  given  fineness  ratio,  the 
shape  calculated  with  the  modified  theory  has  more  curvature  in  the  region 
aft  of  the  nose  than  does  the  body  calculated  with  the  impac.c  theory. 

This  result  is  not  surprising  in  view  of  the  pressure  relieving  effect 
of  centrifugal  forces. 

Calculation  of  the  drag  of  these  bodies  indicates  that  those  obtained 
using  Newtonian  theory  will  as  expected  have  the  higher  drag  at  hypersonic 
speeds,  although  not  by  more  than  a  few  percent.  This  result  suggests 

The  highest  value  of  K  in  the  present  tests  was  2.1. 
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that  consideration  of  centrifugal  forces  will,  in  the  practical  case, 
principally  influence  the  shape  and  not  the  drag  of  minimum  drag  bodies. 


CONCLUDING  REMARKS 


It  has  been  undertaken  in  this  report  to  determine  approximately 
the  shapes  of  several  bodies  having  minimum  pressure  foredrag  at  high 
supersonic  airspeeds.  With  the  aid  of  Newtons  law  of  resistance  and 
the  calculus  of  variations,  an  investigation  was  carried  out  for  various 
combinations  of  the  conditions  of  given  body  length,  base  diameter,  sur¬ 
face  area,  and  volume.  In  general,  it  was  found  that  when  the  length  is 
fixed,  the  body  has  a  blunt  nose  (i.e.,  a  finite  area  of  infinite  slope 
at  the  nose)  as  in  the  classical  problem  considered  by  Newton;  whereas 
when  the  length  is  not  fixed  the  body  has  a  sharp  nose. 

Several  bodies  of  revolution  of  fineness  ratios  3  and  including 
the  calculated  minimum-drag  bodies  for  given  length  and  base  diameter 
and  for  given  base  diameter  and  surface  area,  were  tested  at  Mach  num¬ 
bers  from  2.73  to  6.28  in  the  Ames  10-  by  1 4- inch  supersonic  wind  tunnel. 

A  comparison  of  the  relative  theoretical  and  experimental  foredrag  coef¬ 
ficients  indicated  that  the  calculated  minimum-drag  bodies  were  reasonable 
approximations  to  the  correct  shapes.  It  was  verified,  for  example,  that 
the  minimum-drag  body  for  a  given  length  and  base  diameter  has  as  much  as 
20  percent  less  foredrag  than  a  cone  of  the  same  fineness  ratio.  The  cone 
is,  however,  the  calculated  minimum-drag  body  for  a  given  base  diameter 
and  surface  area. 

The  comparison  between  theory  and  experiment  also  indicated  that  the 
centrifugal  forces  in  the  flow  about  bodies  curved  in  the  stream  direction 
may  influence  their  drag.  The  relative  extent  of  this  influence  was  found 
to  be  predictable,  particularly  at  the  higher  Mach  numbers,  with  a  simple 
modification  to  the  impact  theory  of  Newton.  It  was  therefore  suggested 
that  improved  approximations  to  minimum  foredrag  shapes  at  high  supersonic 
airspeeds  (for  which  the  hypersonic  similarity  parameter  has  a  value 
appreciably  greater  than  l)  may  be  calculated  with  the  aid  of  the  modified 
impact  theory.  Such  a  calculation  was  carried  out  for  bodies  with  the 
same  given  conditions  as  those  calculated  with  the  Newtonian  theory.  In 
general,  the  resulting  shapes  were  found  to  be  somewhat  blunter  in  the 
region  of  the  nose,  to  have  more  curvature  in  the  region  downstream  of 
the  nose,  and  to  have  slightly  lower  drag  than  the  corresponding  shapes 
obtained  using  the  simple  impact  theory. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Dec.  l4,  1955 
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APPENDIX  A 

ESTIMATED  EFFECT  OF  CENTRIFUGAL  FORCES 
ON  SURFACE  PRESSURE  COEFFICIENTS 


An  estimate  of  the  effect  of  centrifugal  forces  on  the  pressures 
at  the  surface  of  a  body  operating  at  high  but  finite  Mach  numbers  may 
be  obtained  by  comparing  the  disturbance  flow  fields  at  these  Mach 
numbers  with  that  associated  with  infinitely  large  Mach  number. 

At  high  Mach  numbers  the  disturbed  air  flows  in  a  relatively  thin 
region  (sometimes  termed  the  hypersonic  boundary  layer)  between  the  bow 
shock  wave  and  the  surface  of  the  body  (see  sketch). 


The  change  in  pressure  from  the  surface  to  the  shock  due  to  centrifugal 
forces  in  the  fluid  is  given  by  the  equation 


assuming  the  directions  of  the  normals  to  the  streamlines  between  the 
surface  and  the  shock  do  not  differ  appreciably  from  the  direction  of 
the  normal  to  the  surface.  This  expression  is  more  conveniently  written 
in  the  form 

-  rN 

AP  -  E  /  P  UdN  (Al) 

RJo 


where  U  and  "E  are  mean  values  of  the  velocity  and  radius ,  respec¬ 
tively,  in  the  interval  N.  Now  the  mass  m  of  air  between  the  surface 
and  the  shock  flowing  (in  unit  time)  by  a  point  on  the  body  is  given  by 
the  relation 


m 


2iry 


L 


N 


~  .  2 


p  UdN  =  jty‘ 


Po  uo 


(A2) 
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Combining  equations  ( A1 )  and  (A2 )  there  is  then  obtained  for  the  pres¬ 
sure  change 


or  in  coefficient  form 


ACp 


££  =  1.2- 

q°  B 


(A3) 


Now  in  the  limit  as  the  Mach  number  approaches  infinity  and  7  of  the 
disturbed  fluid  approaches  1,  the  thickness  of  the  layer  becomes  infini¬ 
tesimal  and  hence 


E  -Eb 


(A*) 


Similarly,  it  is  easily  shown  (e.g. ,  with  the  compatibility  equations 
applying  along  characteristic  lines  in  axially  symmetric  supersonic 
flow)  that 

dTJ  =  0 

along  any  streamline  downstream  of  the  bow  shock,  and  thus  that 

20 

IT  *  — §  /  y  cos  6  dy  (A5) 

y  Jo 


Hence,  in  this  limiting  case,  equation  (A3)  takes  on  a  form  equivalent 
to  that  first  deduced  by  Busemann  (see  second  term  on  right  of  equa¬ 
tion  (2)),  and  later  derived  in  reference  18,  namely, 

rJ 

ACp  *  /  y  cos  6  dy  (A 6) 


where 


sin  6 

dy 


On  the  other  hand  when  the  Mach  number  is  finite,  but  high,  and  7 
of  the  disturbed  fluid  is  closer  to  l.U  than  1,  the  preceding  evalua¬ 
tions  of  E  and  tT  are  in  considerable  error  since  the  hypersonic 
boundary  layer,  although  thin,  is  no  longer  of  infinitesimal  thickness. 
This  change  in  the  boundary  layer  results  from  the  fact  that  the  bow 
shock  is  detached  (except  perhaps  at  the  nose)  from  the  surface  of  the 
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1 


1 


( 

{  , 


\ 


body,  the  lateral  distance  from  the  surface  to  the  shock  increasing 
with  increasing  distance  downstream  from  the  nose  (see  sketch).  Thus, 
for  example,  R  would  he  expected  to  approach  Rg  only  near  the  nose, 
while  with  increasing  distance  downstream  of  the  nose  it  would  he 
expected  to  become  larger  than  Kb*  From  the  pressure  distributions 
presented  in  reference  13  it  is  indicated ,_in  fact,  that  for  K>1  (the 
range  of  K's  of  interest  in  this  paper)  R»Rb  near  the  maximum 
ordinate  of  the  body.  (This  indication  follows  from  the  small  values 
of  the  pressure  coefficients  near  the  maximum  ordinate. )  It  is  sug¬ 
gested,  therefore,  that  at  the  high  supersonic  speeds  under  consider¬ 
ation,  an  approximation  to  R  is  given  by  the  relation 


(AT) 


Similarly,  in  the  case  of  U  it  no  longer  follows  that  the  magni¬ 
tude  of  the  velocity  must  be  constant  along  streamlines  downstream  of 
the  bow  shock  since  pressure  disturbances  can  now  be_transmitted  across 
streamlines.  Thus  a  better  first  approximation  to  U  than  that  given 
by  equation  (A5)  may  be  obtained  from  the  simple  corpuscular  or  impact 
theory,10  namely. 


=  U0  cos  5 


(A8) 


Combining  equations  (A3),  (A7),  and  (A8),  the  estimated  change  in  pres¬ 
sure  coefficient  at  the  surface  of  a  body  due  to  centrifugal  forces  in 
high  supersonic  speed  flow  is  obtained  in  the  form 


or 


ACp  =  ( 1  —  -^-)  cos  5 

% 

ACP  =  £  (1  -  r-)  (sin28)  (A9) 

*2  72  dy 


10With  this  theory,  acceleration  of  the  flow  along  a  streamline  is 
qualitatively  accounted  for. 


i 


,  _  ,  ,  m*  r"  '  *  "  r 
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APPENDIX  B 

CALCULATION  OF  MINIMUM-DRAG  BODIES,  WITH 
CONSIDERATION  OF  CENTRIFUGAL  FORCES 
IN  THE  DISTURBED  FLOW  FIELD 
Given  Length  and  Base  Diameter 

For  the  purpose  of  this  calculation,  equations  (l)  and  (A9)  for 
the  pressure  coefficient  are  combined  with  equation  (4)  to  yield  the 
drag  parameter  In  the  form 


ip  *  yi2  +  <P  (yi)  + 


2  sin2  5  +  2,  (1_  JL)  -jL  sin2  &J  yy'dx  (Bl) 


The  term  yi2  represents  the  drag  on  any  finite  region  of  infinite 
slope  at  the  nose,  vhile  the  function  cp(yi),  given  by11 

<p(yi)  *  -  d  -  cos2  &i 

represents  a  "leading-edge  thrust"  due  to  the  acceleration  of  the  air 
flow  about  a  comer  (if  it  exists)  at  (o,yi). 

The  expression  (Bl)  may  be  put  In  the  form 


ID  "  7iS+ j  |  [2  sin2  S  +  |  (1-  X)  sin2  &J  yy'  -  £  <P(y)  j-  dx 


whereupon  (recalling  sin2  S  =  ■  -  — — g)  the  integrand  simplifies  to  a 

1+y1 

function  f  given  by  the  relation 


11This  function  may  be  obtained  by  evaluating 


T  (1  -  fe  (sin2  5)dy 

yi-€ 

along  the  body-surface  streamline  about  the  comer  at  (o,yi  1 


2b 
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f  =  yy*  2  - 


1  +  lx 

i+2y2 


1  +  y 


»2 


With  the  aid  of  this  expression  and  equations  (15)  and  (19)  the  para¬ 
metric  representation  of  the  minimizing  curve  can  be  obtained  in  the 
following  form: 


y 

y  3 


-1  +  /I  + 


x  = 


i5y2,a  (1  +  y,2)‘ 


(1  +  y2,2)‘ 


,t3 


<iy 

F 


(B2) 


where 


yi*  =  i 


and,  in  general. 


yi  >  0 

The  minimizing  curve  given  by  these  relations,  similar  to  the  curve 
obtained  from  the  impact  pressure  treatment,  does  not  have  a  corner 
between  the  points  (o,  yx)  and  (x2,  y2) .  The  minimum-drag  shape  defined 
by  equation  (B2)  is  compared  in  figure  8  with  that  determined  earlier  by 
considering  impact  pressures  only. 

The  equations  defining  the  minimizing  curves  for  -the  other  given 
geometric  conditions  are  obtained  in  a  similar  manner. 


Given  Length  and  Volume 


-y'3  + 


y  = 


J 


y'6  +  o(l  +  y'2) 


-  Ml  +  y*2)2' 

J2 


-  A(l  +  y,a)‘ 

^2 


(B3) 
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with  y±'  =1  and  y2*  =  0.27^  and  a  range  of  A  from  -00  to  0  correspond¬ 
ing  to  the  range  0  to  *  for  l . 

Numerical  integration  of  equation  (B3)  is  accomplished  by  first 
evaluating  the  first  integral  to  the  Euler  equation  at  the  base  of  the 
body  and  solving  for  c/y2  in  terms  of  y2A.  Letting  cp(y* ,  y2A)  repre¬ 
sent  the  resulting  function  of  y*  and  y2A,  equation  (B3)  becomes 

y  =  y2<p(y*,  ya>0 


'yi 

and  the  volume  is  given  by 


r y2  ,  r O  .274 

i  =  L  #  =  yaA(yaX) 


v  _  v  3 
*  -  y 2 


> 0  .274 


'y*  =i 


q*2  $  =  y23r(y2x) 


The  values  of  the  functions  A  and  T  are  obtained  by  numerical  integra¬ 
tion  for  various  values  of  y2A  to  enable  interpolation  for  that  value 
of  y2A  which  makes  r/A3  =  V/rtZ3.  The  set  (y2A,  A,  r)  so  determined 
satisfies  the  given  volume  and  length  requirements  and  yields  a  value  of 
the  base  ordinate,  y2  =  z/A. 


Given  Length  and  Surface  Area 


y  = 


6y'~ 


[x(i  +  y’2)3'2  -  2y,3]y2  ± 


J  A(l  +  y,2)3/2-2y'3  2y22  +  12cy’3(l  +  y,2)2y2 


>  (b4) 


X  = 


'O 


dy 

y1 


with  a  range  of  A  given  by 

-0.64  <  A  <  2(1  +  ||) 

The  procedure  used  to  integrate  equations  (b4)  is  similar  to  that  employed 
to  integrate  equations  (B3)  above. 
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Given  Base  Diameter  and  Volume 


where 


y  = 


2y2 


Wi  +  y,2)‘ 

•vr*3 


-  3 


j 


*yz  =  5ys,3/(i  +  y2*2)2 


(B5) 


For  i/d  ratios  greater  than  l/2 ,  yx  =  0,  yx*  =  0,  and  the  ranges  of 
y2*  and  A  are 


0  <  y2*  <  «/3 


Given  Base  Diameter  and  Surface  Area 


A  =  121.6(y22/S)3 

The  minimizing  curve  given  by  equations  (B6)  is  compared  in  figure  9  with 
that  determined  earlier  (the  cone)  by  considering  impact  pressures  only. 
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Figure  tr  Comparison  of  approximate  and  exact  pressure  distributions 
over  a  tangent  ogive  of  fineness  ratio  3  operating  at  a  Mach 
number  of  6  (K*2). 
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Figure  2 —  Minimum  drag  bodies  for  various  given  conditions.  (t/ds5.0  for  all  bodies) 


Figure  3-Comparison  of  minimum  drag  bodies  of  given  length  and  base  diameter  determined 

by  tinea:  theory  and  by  impact  theory. 
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tion  for  given  lengths  and  base  -diameters  with  the  approximate 
profiles  employed  in  the  present  tests. 


' 


Figure  5»”  Photographs  of  the  seven  test  bodies,  the  shapes  of  which  are  given  by  the  equation 

y  =  (x/z)n,  and  the  if d  =  3  tangent  ogive. 
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(b)  / Vd  *  5  bodies 

Figure  7-  The  ratios  of  foredrag  coefficients  of  test  bodies  to  foredrag 
coefficients  of  cones  as  functions  of  the  exponent,  n,  in  the  equation 
defining  body  shapes. 


Figure  9-  The  effect  of  centrifugal  forces  on  the  shape  of  the  minimum  drag  body 
of  given  diameter  and  surface  area  (d  =  2,  5  s  31.57). 
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